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The intermediate coupling regime in polaronic systems, situated between the adiabatic and the 
anti-adiabatic limit, is characterized by resonant pairing between quasi- free electrons which is in- 
duced by an exchange interaction with localized bipolarons. The onset of this resonant pairing takes 
place below a characteristic temperature T* and is manifest in the opening of a pseudogap in the 
density of states of the electrons. The variation of T* is examined here as a function of (i) the 
typical frequency ujo of the local lattice modes, which determines the binding energy of the bipo- 
larons, and (ii) the total concentration of charge carriers ntot = np +nB, where np and ub are the 
densities of free electrons and bipolarons, respectively. The variation of either of these quantities 
induces similar changes of the value of ns with respect to that of wf, in this way leading to a shift 
of r*. For finite, but small values of ub (< 0.1 per site), we find negative and practically doping 
independent values of the corresponding isotope coefficient a*. Upon decreasing ntot such that ub 
becomes exponentially small, we find a rapid change in sign of a*. This is related to the fact that 
the system approaches a state which is more BCS-like, where electron pairing occurs via virtual 
excitations into bipolaronic states and where T* coincides with the onset of superconductivity. 

PACS numbers: PACS numbers: 74.20.Mn, 74.25.-q, 74.25. Kc 



I. INTRODUCTION 

The experimental verification of an isotope effect in 
the classical low temperature superconductors has been 
an unequivocal proof for phonon mediated electron pair- 
ing in form of Cooper pairs. Such a clearcut proof is 
so far absent in the high temperature superconducting 
cuprates (HTSC). Moreover, their proximity in the un- 
derdoped regime to an antiferromagnetic insulating state 
has largely contributed to conclude that in those mate- 
rials pairing is associated with strong correlations. Nev- 
ertheless, the question concerning the origin of the pair- 
ing is far from being settled and lattice driven pairing 
ought not to be ruled out at this stage of our understand- 
ing. Actually, there is a certain amount of experimen- 
tal evidence that strong electron-lattice coupling plays 
some role in stabilizing the superconducting phase in the 
cuprates. What is clear is that the pairing is definitely 
not of the form of a phonon mediated BCS one, corre- 
sponding to a weak coupling adiabatic regime. 

Early on, anomalous mid-infrared optical absorption 
was found in virtually everyone of the superconducting 
cuprates. In lanthanum based compounds the supercon- 
ducting critical temperature Tc was found to scale with 
the oscillator strength of this absorptioni*^ and later on 
was shown to be due to polaronic charge carriers^. More 
refined measurements followed. From neutron scatter- 
ing studies it became clear that the manifestation of 
the superconducting state in the cuprates could be in- 
timately linked to strong electron-lattice coupling. A de- 
tailed analysis of the phonon density of states showed 
that the high-frequency modes are significantly renor- 



malized in the superconducting materials as compared 
to their insulating parent compounds^. A further man- 
ifestation of strong electron-lattice coupling comes from 
the observation of a kink in the electron quasi-particle 
dispersion in the 50-80 mev energy region seen in angle 
resolved photo emission spectroscop}^. Finally, inelastic 
neutron scattering experiments pointed to an anomalous 
behavior in the dispersion of the in-plane longitudinal op- 
tical phonons with wave vectors [0, 0.25, 0] in the YBCO 
superconductors. This corresponds to bond stretching 
vibrations being associated with dynamical charge fluc- 
tuations on the Cu ions driven by the displacement of 
the neighboring ligand environment of the O atoms^ . 

Let us suppose, as a working hypothesis, that su- 
perconductivity in the cuprates is indeed controlled by 
strong electron-lattice coupling. If we want to test this 
assumption by examining the isotope effect in those ma- 
terials, the right quantity to look at is not the transition 
temperature Tc, but the onset temperature of the elec- 
tron pairing, T*. T* shows up in a qualitative change 
of the photoemission spectrum such that the electronic 
density of states exhibits a charge pseudogap as the tem- 
perature is decreased below T* , eventually merging into 
a true superconducting gap below Tc- Pair formation in 
BCS superconductors coincides with the onset of a global 
phase-coherent superfluid state and hence the isotope ef- 
fect can be evaluated on the basis of the shift in Tc- This 
does not apply to the HTSC, where it is a pair resonance 
state which sets in below T*, implying pairing on a fi- 
nite length and time scale. Only when this length and 
time scale gets longer and longer upon decreasing the 
temperature, a global phase-coherent state can be estab- 
lished, which is controlled by the center of mass motion 
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of the Cooper pairs rather than by their breaking up into 
individual electron pairs, as in the case of BCS supercon- 
ductors. Considering that electron pairing is of resonance 
type rather than of a true bound electron pair nature, the 
different experimental setups devised to capture such a 
feature must rely on a time scale short enough to see 
this pairing as static. Thus, NMR or NQR cannot de- 
tect it since the relevant time scale has to be well below 
10^^ sec. On the contrary neutron spectroscopy, study- 
ing the relaxation rate of the crystal field excitations, 
and X-ray absorption near edge spectroscopy (XANES) 
are in the right time scale regime of [10~^^, 10""'^^]. And 
in fact it is those measurements on La2_a;Sra;Cu042*i£, 
HoBa2Cu40s^ and Lai.8iHoo.o4Sro.i5Cu04^ which ini- 
tially demonstrated this resonant pairing isotope effect. 

We shall in this paper explore the resonant pairing 
isotope effect on the basis of a phenomenological model, 
believed to capture the intermediate coupling regime in 
polaronic systems, situated between the standard weak 
coupling adiabatic Born-Oppenheimer regime (applicable 
to BCS superconductors) and the anti-adiabatic regime, 
where the electrons pair up into bipolarons, expected to 
be localized at low temperatures. We have previously in- 
troduced such a model to study the connection between 
local dynamical lattice deformations, measurable by EX- 
AFS pair distribution functiong^^, and the incoherent 
background expected in photoemission spectra, as inter- 
preted in terms of phonon shake-off processesi^. On the 
basis of this model we shall show that the isotope effect 
of T* can be traced back to the pairing energy of the 
bipolarons which is a linear function of the characteristic 
frequency of the local lattice modes. 

In section II we briefly sketch this model and present 
a scheme to determine the doping and frequency depen- 
dence of T*. Section III is devoted to a discussion of the 
results on the isotope coefficient and an attempt is made 
to relate its behavior to the isotope shift of the bipo- 
laron binding energy and to the renormalization of the 
exchange coupling between the localized bipolaron and 
itinerant electrons. In the Conclusions, section IV, we 
make a comparison of these results with what is known 
presently from the the experiments and suggest new ways 
of looking at this problem. 



II. THE MODEL AND RESONATING PAIRING 
TEMPERATURE 

Exact diagonalization studies^^ in the crossover regime 
between the adiabatic weak coupling limit and the anti- 
adiabatic limit have led us to the conjecture that in this 
regime we are facing strong fluctuations between tightly 
bound pairs and uncorrelated pairs of free electrons, 
a scenario which can be phrased into a boson-fermion 
model. In order to incorporate into this model the infor- 
mation concerning the origin of the tightly bound elec- 
tron pairs, we assume explicitly that they are of bipola- 
ronic nature. The minimum model which can describe 



such a situation is then given by the Hamiltonian 



H = {D - fl)J2n^a -t cf^Cja 
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Here pf denote the creation and annihilation operators 
for the electron pairs which, due to their interaction with 
the local lattice deformations Xi = {ai + al)/ \j2MiJo/Ti^ 
end up in self-trapped bipolarons pf exp [±Q!(ai — af)\ 
localized on some effective sites i. Such entities are 
treated as hard-core bosons with spin-i commutation re- 
lations, = 2pf and [pt,P'j]^ = a^t^ de- 
note annihilation (creation) operators of the excitations 
of local lattice displacements, M is some atomic mass 
characterizing the effective sites, and wq is the frequency 
of the dynamical local lattice deformations, cf'^ are the 
annihilation (creation) operators for the itinerant elec- 
trons with spin cr, rijo- = cf^Ci„ being the related number 
operator. The bare nearest-neighbor hopping integral for 
such electrons is given by i, corresponding to a band- 
width 2D = 2zt where z denotes the lattice coordination 
number. The other parameters of the model are the bare 
bosonic energy level A b , the coupling a of the local elec- 
tron pairs to the surrounding lattice deformations and 
the bare exchange coupling v between the bosons and 
the pairs of itinerant electrons. The chemical potential 
/i, being common to electrons and bosons, guarantees the 
overall charge conservation. This model has been studied 
extensively in the limit of zero coupling (a — 0) to the 
lattice, particularly in connection with the pseudogap ef- 
fect in the HTSC induced by the resonant pairing of the 
electrons. The basic idea behind it is that the closeness of 
a weakly bound two-electron state to the energy level of 
two itinerant electrons induces resonant pairing with sub- 
stantial lifetime in the electronic subsystem. This is an 
effect which is analogous to the atom pairing induced by 
Feshbach resonance in trapped ultracold gases, studied 
in connection with their condensed statea^^. The onset 
temperature T* for electron pairing is determined by the 
strong drop-off with decreasing temperature of the on- 
site correlation function {pfciiCi]) when passing through 
T — T*—. The change in this correlation function is in- 
dependent on any onset of long-range phase coherence 
and is described by mere amplitude fluctuations. Using 
a variational wave function of the form 
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FIG. 1: T* as a function of ritot for a variety of different local 
phonon frequencies uoq 



describes the pseudogap regime in an approximate form. 
For the purpose of the present study we shall adopt such a 
mean field analysis for which we shall here merely sketch 
the procedure. 

Since the mean-field decoupling leading to Eq. (3) sepa- 
rates the fermionic part from the bosonic bound electron 
pairs, we can write the eigenstates of this Hamiltonian 
as a direct product of the two separate Hilbert spaces 
associated with fermions and bosons in the form 
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with 



l^^> = T{[uk+vkcl^ct^^] |0) 
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(7) 



we see that the exchange coupling term of the Hamilto- 
nian (1) becomes 



vpx -f vx 



with 
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denoting the amplitudes of the order parameters of the 
electron and boson subsystems. Thus, any inter-site 
phase fluctuations are explicitly suppressed in such a 
mean-field approximation. It hence guarantees that the 
resulting transition is exclusively due to amplitude fluc- 
tuations and thereby lends itself to describe the onset of 
pairing without any simultaneous onset of phase coher- 
ence. That this approach for determining T* is qualita- 
tively and, to a large extent, also quantitatively correct, 
was checked with a comparative study based on exact 
diagonalization proceduresi^ and self-consistent pertur- 
bative approaches'^. The rapid but smooth drop-off of 
the local correlation function (p^CiiCi^) at T* in those 
studies is then apparent as a sharp drop-off to zero of 
the same function at a mean-fleld critical temperature 
"^MFA- Generally T^pj^ is found to lie slightly below 
T* determined in more elaborate treatments, but shows 
the same dependence on the charge carrier concentra- 
tion. This justifies the use of an analogous mean-field 
type procedure for the generalized boson-fermion model 
presented above in Eqs. (1-4), in order to extract a T* 
when the coupling of the lattice vibrations to the charge 
carriers is turned on. 

A detailed account of this mean-field analysis has been 
given in ref. [11] where we associated our results to a 
superconducting phase, assumed to be controlled exclu- 
sively by amplitude fluctuations. The true superconduct- 
ing phase for this model is however known to be con- 
trolled by phase fluctuations'^, while its mean-field phase 



Here |0) and |0) are the vacuum states for fermions and 
bosons, respectively, and \n) is the n-th excited harmonic 
oscillator state. It should be noted that phonons are only 
connected to bosons and thus the contribution (7) is the 
only one requiring numerical diagonalization. Denoting 
the eigenvalues of the two mean-field states (6) and (7) 
by £k(p) ~ ±\/(ek — A*)^ + (^'p)^/4 (which differs from 
the bare electron dispersion by showing a gap of size 
vp) and Ei{x), we have the following selfconsistent equa- 
tions for the order parameters and the concentration of 
electrons and local pairs 



1 tanh^ 

^ UlnVln exp [-PEi{x)\ 



(8) 
(9) 



ntot 



tanh^ 



(vin)^] exp i-lSEii^)]- (10) 



Here Z = exp [— /3i?;(a;)] denotes the partition fmic- 
tion corresponding to the bosonic part of the mean-field 
Hamiltonian, given by 

Hb = (Ab - 2m) J2 (p^ ^^T.^Pt + A"] 

-hujoa ^ (p^^ + i j (a, + af) + fiujo ^ afa., . (11) 

i ^ ^ i 

The onset temperature T* for electron pairing is then 
determined by solving the above set of equations in the 
limit a: — > 0, p — > 0. 

We assume an electronic band extending from —D to 
D and choose a set of parameters a = 2, = 0.1, v = 
0.25, with the phonon frequency ujq varying in the range 
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ntot 

FIG. 2: T* (top panel) and 2nB (bottom pannel) as a function 
of ntot for a variety of different local phonon frequencies loq 



[0.01,0.1], such as to cover the intermediate polaronic 
situation (all energies arc in units of the half-bandwidth 
D). This choice, leading to values of T* of the order of a 
few hundred degrees K (T* ~ lO^D), ensures that upon 
changing the total number of charge carriers ntot, one 
covers the regime of electron concentration close to half 
filling, with the possibility of having a drastic decrease of 
T* with small variations of ntot- 

We present in Fig. 1 the variation of T* as a function 
of ntot for several values of the local phonon frequency 
ujQ. As Wo is increased, we observe the following two main 
effects: (i) a shift of the whole curve T*{ntot) to lower 
values of ntot, and {ii) an overall diminution of the value 
of T* . The first effect is due to a shift (5A b — — £bp 
of the bosonic energy level associated with the bipolaron 
binding energy Ebp — oP'Tiujq. The second effect is due 
to the decrease of the effective exchange coupling term, 
determined by the reduced overlap of the lattice defor- 
mations corresponding to the presence of bipolarons and 
free electrons, respectively, on a given site. In the ex- 
treme strong coupling limit o^Tiloq > D this renormal- 
ization would correspond to w — > ve~°' , but remains of 
reasonable size and is frequency dependent in the inter- 
mediate coupling case. In conclusion, reducing ntot with 
doping or increasing ujq upon increasing the isotope mass, 
leads to qualitatively similar results in the shift of T* as 
already pointed out by some experimental observations^. 



III. THE RESONATING PAIRING ISOTOPE 
EFFECT 

We shall focus here on a regime of doping (ntot) and lo- 
cal phonon frequency (wq) where T* shows a rapid drop- 
off upon decreasing the total number of charge carriers 




FIG. 3: T* (top panel) and 2nB (bottom panel) as a function 
of the local phonon frequency ljq and for a variety of different 

ntot 



ntot over a relatively small range of values, i.e., [0.8, 1.2]. 
This choice is made in an attempt to account for the 
anomalous behavior of the isotope coefficient a* observed 
in the HTSC, which in the underdoped to optimal dop- 
ing regime shows unusual and negative values. In order 
to show the evolution of this isotope effect upon going 
from the underdoped to the overdoped case (the latter as- 
sumed to be more BCS-like), we examine how a* changes 
with the particle density. For doping rates ntot such that 
ns is finite but small we have resonant pairing of the 
electron pairs due to their exchange with localized bipo- 
larons. On the contrary, if is exponentially small these 
exchange processes are purely virtual processes. We then 
have a situation where two-particle pairing can no longer 
be realized except via a true many-body effect describ- 
ing Cooper pairing. The corresponding T* then signals 
simultaneous pairing and onset of superconductivity to- 
gether with an isotope exponent which is positive. 

The behavior of T* as a function of ntot and ujq is 
controlled by two competing effects: 

(z) T* increases with increasing ntot when ns varies 
between and 0.5. 

(ii) T*, for a fixed ns, decreases with increasing ujq be- 
cause of the polaron induced reduction of the exchange 
coupling V. But since an increase of ujq not only reduces 
the effective exchange coupling but also leads to an in- 
creased bipolaron level shift which in turn increases ns 
for a fixed ntot, T* is in general a non-monotonic function 

of LUq. 

We present in Fig. 2 the variation of T* and ns as a 
function of ntot in the above mentioned regime of pa- 
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rameters. We notice the onset of a rapid rise of T* with 
increasing ritot controlled by an equally rapid rise oing. 
Upon further increasing ntot^ T* starts to saturate, an 
effect due to the hard-core nature of the localized bipo- 
larons which becomes important when ns approaches 
0.5. 

The variation of T* and ub with ljq for a set of different 
ritot is illustrated in Fig. 3. For sufficiently large values 
of ritot, such that ub is finite but still lower than 0.5, we 
find a monotonically decreasing function with increasing 
t^Oj controlled by the polaron induced reduction of the 
effective exchange coupling v. As ntot decreases, leading 
to a vanishing concentration of bipolarons, the behavior 
of T* changes qualitatively. For low values of cjo, the 
polaron induced reduction of the effective exchange cou- 
pling is negligible and T* now increases with increasing 
luq because it is controlled by the increase of ub- With 
increasing wq, the effect of the polaron induced reduction 
of the exchange coupling becomes competitive with the 
increase of ub such that the initial increase of T* with 
increasing wg changes into a decreasing behavior. 

Provided we are in the regime of resonant pairing, with 
Ub small but finite and T* monotonically decreasing with 
decreasing ntot, Fig- 3 tells us the following. T* is shifted 
upwards by the decrease of the phonon frequency asso- 
ciated with the increase of the isotope mass (as realized, 
for instance, replacing ^^O by ^'^O), and this effect be- 
comes less and less pronounced as doping increases (ritot 
decreases). Moreover, if for a given doping level the iso- 
tope substitution is made for heavier elements, such as 
^^Cu replaced by ^^Cu, the increase of the correspond- 
ing T* is getting smaller and smaller, as one can deduce 
from the behavior of T* at low ojq. These features are 
in qualitative agreement with experimental findings in 
LaSrHoCu04 compounds^iiSi. 

In order to determine the isotope coefficient, which it- 
self will depend on loq and ritot, we have first interpolated 
the calculated values of T* by a ratio of polynomials, and 
then derived a* from the relation 

a = 0.5 In / — — — - (12 

\T*{i + l) cJo(« + l)/ 

where wo(«) and wo(j -I- 1) represent two very close values 
of the phonon frequency. We present in Fig. 0] the varia- 
tion of the isotope exponent as a function of the phonon 
frequency for a set of different values of ritot and in 
Fig-El the variation of the isotope exponent as a function 
of ritot for a selected set of phonon frequencies ujq. 
We notice two distinct regimes which characterize a*: 

(i) a regime where a* is negative and depends strongly 
on LjQ but shows a relative independence on the concen- 
tration ritot- This happens when the corresponding is 
big enough to sustain the mechanism of resonant pairing. 

(ii) a regime of positive values of a* which occurs when 
n B drops to zero as a consequence of the decrease of ritot 
or, alternatively, of wq (see the bottom panel of Fig. 2). 
In this case, which corresponds to the system being more 
BCS-like, electron pairing arises from virtual excitations 
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FIG. 4: Q* as a function of the local phonon frequency cjq 
and for a variety of different ntot 
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FIG. 5: Q* as a function of ntot and for a variety of different 
local phonon frequencies luq- 

of the electrons into bipolaronic localized states above 
the Fermi level and the onset of pairing at T* coincides 
with the onset of superconductivity. 

These two regimes are clearly visible in the variation 
of the gap ratio 2A{0)/{kBT*), where 2A(0) = vp{0) 
is the zero temperature gap in the Fermionic excitation 
spectrum ik- This ratio, which is shown in Fig. 1^1 as a 
function oiritot for several ujq, a, slowly- varying function 
of the concentration when the system is well inside the 
resonant pairing regime (1.2 < ntot < 2.8 for our choice 
of luq values) but strongly depends on the local phonon 
frequency loq. For values of ntot such that becomes 
exponentially small (depending on this happens below 
ntot — 0.9, see bottom panel in Fig. 2), the gap ratio 
approaches the BCS value 3.52. 

Going back to the behavior of the isotope coefficient 
at low ntot, we stress that even in this BCS-like regime 
a* shows a frequency dependence related to the pair- 
ing mechanism, which in the case studied here of strong 
electron-phonon interaction is different from the standard 
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FIG. 6: 2A(0)/(fc_Br*) as a function oi ntot and for a variety 
of different local phonon frequencies ujq. 



weak-coupling BCS one. In our model, pairing among 
the electrons is induced via virtual excitations of elec- 
tron pairs into a localized bipolaron state and an effective 
BCS-like Hamiltonian for this situation can be derived 
along the following two steps: 

i) the boson-phonon coupling is incorporated into 
an effective boson-fermion exchange interaction v — 
veicp (— via a Lang-Firsov approximation; 

ii) the boson-fermion coupling term is subsequently 
eliminated to linear order via the usual unitary trans- 
formation 

H = e^iJe"^ (13) 
S = ^/(fc)[p+c_fe^CfeT-H.c.] (14) 



i,k 



fik) 



2£k - As + EBP 



(15) 



which in the end results in an effective BCS-like Hamil- 
tonian of the form 



Y E [f{k) + f{k')]cl^c+^^c.u'icu'^ .(16) 



k,k' 



In the standard weak-coupling effective BCS Hamilto- 
nian, the interaction term is restricted to the summation 
over fc-vectors within a small energy range around the 
Fermi surface of width given by the Debye temperature 
Qd. This gives rise to an effective interaction of the or- 
der of from which a frequency independent isotope 
shift follows, with an isotope exponent equal to 0.5. 

The effective Hamiltonian for our scenario, eg . 1161 
shows on the contrary no such a cutoff in energy and 
extends the pairing of electrons over all k vectors below 
the Fermi vector, albeit with different weight. Moreover, 
because of the presence of the bipolaronic energy in the 



denominator of eg . 1151 the effective interaction is attrac- 
tive in this regime of low bipolaron concentration, since 
2ekp < As ~ Sbp, and grows in magnitude as ojq in- 
creases. This is the reason why we obtain a frequency 
dependent isotope shift even in this BCS-like regime with 
a T* which increases as ujq increases. 



IV. CONCLUSIONS 

Experimentally, there are strong indications that in the 
HTSC the isotope coefficient associated with the tem- 
perature T* at which the pseudogap in the underdoped 
regime opens up, is negative. Its precise numerical value 
depends on the type of material, the doping regime, the 
type of isotope substitution (^^O ^ ^^O , ^^Cu ^ ^^Cu) 
as well as on the different time scales of the various ex- 
periments. Presently, there exist to our knowledge no 
systematic experimental studies which would permit to 
test a particular theoretical approach on this issue in 
detail. The study presented here was designed to in- 
cite experimental work to explore specifically the doping 
dependence of the resonant pairing isotope effect, given 
the characteristic strong doping dependence of T* in the 
HTSC. If we assume that this feature is related to a two- 
component scenario with localized charge carriers (such 
as bipolarons) and itinerant electrons, with the sharp 
drop of T* being essentially governed by a doping de- 
pendent change of the bipolaron thermal population, the 
isotope effect on T* should exhibit a corresponding con- 
centration dependence. In agreement with experiments, 
we find indeed a negative value for the isotope exponent, 
as long as pairing is assured by a resonant scattering 
between the localized bipolarons and the free electrons. 
However, as soon as upon doping we enter a regime where 
the bipolaron level moves above the Fermi energy, pair- 
ing is only possible via a collective effect such as Cooper 
pairing and the isotope exponent switches sign and be- 
comes positive. With respect to the HTSC, this could 
happen when going from the underdoped into the over- 
doped regime. Considering, however, that the lattice me- 
diated coupling between electrons in the present scenario 
is different from the standard phonon-induced Cooper 
pairing, we obtain a behavior for the isotope exponent 
which deviates from that of standard low temperature 
BCS superconductors. Although it converges to a value 
independent on doping as ritot is reduced such that ns 
becomes exponentially small (see the curves at low luq in 
Fig. 5), it nevertheless continues to sensitively depend on 
the characteristic local phonon frequency ujq. 

The present mean-field type study is expected to qual- 
itatively correctly describe the doping dependence of a* 
in the resonant pairing regime where T* is controlled by 
amplitude fluctuations. This mean-field scheme treats 
the local dynamical atomic displacements as correlated 
to the local density fluctuations between the electrons 
and the bipolarons, which in our model, are responsi- 
ble for the opening of the pseudogap. It is this mecha- 
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nism which is at the origin of the blockage of the crystal 
field excitations which below T* couple to the free charge 
carriers, as observed experimentally in neutron spectro- 
scopic measurements testing the transitions between dif- 
ferent crystal field levels^-?. A more quantitative analysis 
than the mean-field procedure presented here would be 
required in order to account for the dynamical nature of 
the onset of pairing, as seen in experiments with differ- 
ent time scalesii which lead to different absolute values of 
T* but to qualitatively similar doping dependence of T* . 
Such a highly non-trivial undertaking is however beyond 
the scope and purpose of the present work. 

We want to stress that the scenario described here is 
based on a mechanism originating from purely local dy- 
namical lattice instabilities for which experimental ev- 
idence has been accumulating over the past few years. 
Inelastic neutron scattering measurements have shown 
strong compositional dependence of certain optical (half- 



breathing zone-edge) phonon modes, which were linked 
to spatial local charge inhomogeneitiesi^ of small clusters 
and suggest that the lattice is strongly involved in the 
charge dynamics. EXAFS studies tracked such local 
charge inhomogeneities in form of a significant deviation 
from a systematic Pauling-type shift of the planar Cu-0 
bonds, when the bonding mechanism changes from ionic 
to more covalent nature as doping is increased. Similarly, 
very recent tunneling spectroscopic studies^ indicate the 
existence of local charge density modulations involving 
local spatial correlations of four Cu02 unit cells, referred 
to as "squared checkerboard" structures. All these find- 
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